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by 

Rudolf A. Hanel and Frank Bartko 
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SUMMARY 



The interpretation of the data obtained for Venus in the 8-12/x 
atmospheric window clearly requires some knowledge of the 
temperature distribution in the planet's upper atmosphere. For 
this purpose temperatures, net fluxes, and specific intensities 
were calculated for a C0 2 -N 2 atmosphere in radiative equilibrium. 
The wavelength, temperature, and pressure dependence of the C0 2 
absorption, including the excited bands, was considered. Solar 
radiation absorbed in the clear atmosphere influences the tem- 
perature profile distinctly. The effects of various C0 2 concen- 
trations, cloud top pressure levels, and cloud reflectivities on the 
equilibrium temperatures are also discussed. 
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RADIATIVE EQUILIBRIUM IN PLANETARY ATMOSPHERES. 
I. APPLICATION OF THE STRONG LINE ABSORPTION 
LAW TO THE ATMOSPHERE OF VENUS* 

by 

Rudolf A. Hanel and Frank Bartko 
Goddard Space Flight Center 


INTRODUCTION 

The structure of a planet's atmosphere is governed primarily by convective and radiative 
energy transport. On the earth, convective processes preponderate in the troposphere and radi- 
ative processes in the stratosphere. Similarly, the atmospheres of Mars and Venus should exhibit 
zones of convective and radiative transfer. Convection probably characterizes the Venus atmos- 
phere below the cloud level, especially if the clouds consist of dust particles, whereas the region 
above the cloud layer is probably dominated by radiative transfer. The interpretation of infrared 
measurements of the atmospheres of the terrestrial planets must, therefore, be based on the 
theory of radiative transfer and appropriately modified for convection where necessary. The ob- 
jectives of the program discussed in this paper were to establish an analytical tool for the study of 
planetary atmospheres based on the radiative transfer theory. An initial version of this program 
has been applied to the atmosphere of Venus. This application must be used with relatively inaccu- 
rate estimates of physical parameters which influence the solution of the transfer equation. How- 
ever, the program permits a parametric study where the number of possible solutions is ultimately 
restricted by observational constraints. 

The equilibrium temperature distribution of an N 2 -C0 2 atmosphere was calculated for several 
volume concentrations of C0 2 . The wavelength, temperature, and pressure dependence of the 
absorption coefficients was taken into account. The effect of direct and diffusely reflected solar 
radiation absorbed by the atmosphere was also included. An effective surface emissivity was 
introduced as an additional parameter and permitted a deviation from the usual assumption that the 
surface radiates like a blackbody in the infrared. In addition to the radiative equilibrium tempera- 
tures, the specific intensities were calculated. Net fluxes were obtained by angular integration, 
taking into account fully the dependence of the intensity on direction. The simplifying procedure of 
converting parallel beam radiation to an equivalent diffuse radiation field was avoided (Refer- 
ences 1-3). 


‘This report was presented at the Conference on Atmospheric Radiation, Eppley Laboratories, Newport, Rhode Island, and will be 
published in an American Meteorological Society Monograph. 
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The net flux of a planet can be determined independently from the heat budget by considering 
the observed albedo and the appropriate value of the solar constant. This restricts the number of 
theoretically possible solutions, since only those which yield total flux values consistent with the 
albedo can be admitted. The number of solutions can be restricted further by a comparison of 
calculated and measured values of the specific intensity at several wavelength intervals and zenith 
angles. 

In the discussion below, the assumptions and limitations of the calculations are stated and the 
mathematical formulation is presented. An illustration of the effects of the various parameters on 
the temperature profiles follows, but further conclusions will be discussed elsewhere. 


SOLUTION OF THE RADIATIVE TRANSFER EQUATION 

In order to calculate the equilibrium temperatures in a finite, nongrey, optically thick atmos- 
phere, the equation of radiative transfer must be solved, subject to appropriate boundary conditions 
(Reference 4, for example). The assumption of radiative equilibrium holds as long as the atmos- 
phere is stable against convection (Reference 5, for example). Stability, in turn, is maintained as 
long as the temperature gradient is everywhere less than the adiabatic gradient. 

For the particular application considered here, the surface is taken to be the cloud tops and is 
assumed to be the demarcation between regions of energy transport dominated by radiation and 
convection. It is realized that the assumption of a solid surface for the clouds is a severe simpli- 
fication. A more rigorous treatment would apply the theory of radiative transfer for a scattering 
medium (Reference 6, for example). 

For an atmosphere in local thermodynamic equilibrium the equation of transfer is 


dl 


vjf 1 ' T »v) 


dr 


V,) - B V (T) 


( 1 ) 


where B v (T) is the Planck function and i^/x, n v \ the specific intensity of the radiation field. For 
radiative transfer calculations dealing with strong molecular band absorption, the optical thickness 


t can be expressed 




r v{ l v< U *) 


( 2 ) 


where i v is a generalized absorption coefficient and u*/m a pressure- and temperature-reduced 
path length. It is convenient to use the optical path in the vertical direction, t v , and the direction 
cosine instead of the optical path, r fiu) in an arbitrary direction. The appearance of jj. under 
the exponent n has a significant influence on the transfer problem. The value of n is unity for 
weak line or grey absorption and one-half or less for the square root or strong line absorption law 
(Reference 7). In contrast to weak line or grey absorption, the strong line absorption law reduces 
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the thermal coupling between adjacent atmospheric layers and permits a more effective exchange 
of radiative energy over greater optical distances. 

Very general boundary conditions have been selected. At some optical depth, t vc , the atmos- 
phere is assumed to be bounded by a perfectly diffuse reflector of effective emissivity e Vc , related 
to the surface reflectivity by 


E »C “ 1 T VC ■ ( 3 ) 

For positive values of /x this assumption requires the solution of Equation 1 at the boundary t vc 
to be 


I + (t = T ) 
V \ V vc/ 


e vc M T c) + r „c F v'( T Vc) + ~7T~ Z 


VC 


(4) 


The first term on the right represents the thermal emission by the surface of temperature T c . 
The second term arises from the nonblack characteristic of the surface; the downward flux inci- 
dent on the surface is partially reflected back into the atmosphere. The third term describes the 
reflected solar component. 

At the top of the atmosphere where t v = o only solar radiation is assumed incident. For 
thermal radiation outer space is considered a perfect sink. Hence the upper boundary condition is 


K(z< r v = °) 


(5) 


The solution of the transfer equation subject to the boundary conditions and the absorption law 
for t may be verified to be 


and 


T - T 
VC . V 


v( * T v) = £ vc B „( T C )‘ 


B„(t) e 


* ~ T v 

M n 


T - T 
VC V 


+ r e 

VC 


S - vc — Vc 

B„(t) E 2 (t vc ~ t)dt + r vc i e ^ e ^ 


( 6 ) 


K (* M 


, n dt 

B v (t) e " — + — i e 


v r n 

— r * S 


(?) 
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In Equation 6 ij (/a r v ) is the upward intensity whose respective components arise from the 
surface emission, the atmospheric emission by all layers at optical depths greater than r v , and the 
downward fluxes diffusely reflected upwards, ij (/!, r v ) is the downward intensity consisting of the 
atmospheric emission arising from all layers with optical depths less than r v and the incident 
solar radiation. 

The temperature distribution is obtained from a solution of the equation of radiative equilibrium 
which equates the total emitted and absorbed radiation for each volume element. For a stratified 
atmosphere each layer must satisfy the energy balance condition 

m l v e M dV do dv = 
v 

The exponential term represents the attenuation losses due to self-absorption in the layer of finite 
optical thickness. Integration over all directions is required since l v as well as e~ AT/fJ - n is a 
function of Furthermore, since absorption and emission may occur in different parts of the 
spectrum the equality holds only if the integration is taken over all frequencies. This general form 
of the condition of radiative equilibrium applies to any volume element as well as to different forms 
of absorption laws. If At is very small in all directions, the familiar form (References 8 and 9) of 
the radiative equilibrium equation results: 


n 

•A,, J . J\r 


I„ dV do) dv 


(8) 


\ - K Jv Av ■ (9) 

Substitution of the intensities obtained (Equations 6 and 7) into the equation of radiative equilibrium 
leads to a Milne-type integral equation in b v (t) , whose solution is required in obtaining the tem- 
perature distribution. The solution is obtained by iteration beginning from an assumed tempera- 
ture distribution. The net flux is computed after each iteration. When the distribution of the net 
flux is constant within specified limits, radiative equilibrium is achieved and the calculation pro- 
ceeds to the final step, the computation of limb darkening. 

The components of the upward and downward intensities, the net flux, and the condition of radi- 
ative equilibrium will be examined in more detail in the manner required for computer solution. 
Hence a brief description of the coordinate system and the calculation of r v follows. 


COORDINATE SYSTEM 

A plane parallel atmosphere is assumed and the computations are carried out in a 33 layer 
coordinate system (Figure 1). The interface between layer j - l and j is called interface j and 
the quantities u* and r. are computed from the top of the atmosphere (interface one, zero pres- 
sure) to interface j and the corresponding pressure P 3 . In all computations the pressure at 
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interface j = 2 was taken equal to 1CT 3 atm. 
Layer one above this level is included for 
the sake of completeness. The computed 
temperature is only an average value and it 
should not be interpreted that layer one is 
isothermal. For this region of the atmos- 
phere, local thermodynamic equilibrium may 
not exist (Reference 10, for example) and the 
absorption law adopted may not strictly ap- 
ply (Reference 11). In addition, photochemi- 
cal and other atmospheric processes will 
have a strong effect on the temperature in 
this region. For layers between j = 2 and 28, 
pressures are calculated according to the 
recursion formula 


INTERFACE 



32 

33 

•//7//rrr?rm 34 

CLOUD SURFACE 


Pj — 0 (u* and T l - 0) 

P 2 = '0" 3 a«m 



P 

c 


Figure 1 —Coordinate system. 



(10) 


This mode of pressure slicing yields layers of approximately the same geometrical height but of 
increasing optical thickness. To avoid a very thick layer adjacent to the cloud boundary, a linear 
pressure slicing was applied between j = 28 and the surface (j = 34). 


ABSORPTION COEFFICIENTS 

The calculation of the optical thickness requires a value for the absorption coefficient, which 
varies with pressure and temperature, and a value for the mass of the absorbing gas. For pres- 
sures greater than 10 " 3 atm the individual line shape is determined primarily by collision broaden- 
ing. The effect of this factor on atmospheric structure has been studied by several authors 
(Reference 12, for example). For band absorption the concept of optical thickness defined for a 
single frequency is still very useful, but r becomes, in general, a complicated function of u*. For 
all cases studied in this paper, an analytical expression can be derived for this function, based on 
the so-called strong line absorption law (Reference 7). This approximation applies to C0 2 within 
the pressure ranges of interest, namely from 10' 3 to 1 atm for C0 2 concentrations exceeding 0.01 
part per volume. 

The mass of active gas is usually expressed in terms of a path length. An element of the path 
length is defined by the product of the geometric path length ds and the normalized density of the 
active gas 


du 



(ii) 
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Since p = MP/RT , 


du q -rf ds . (12) 

In the domain of collision broadening, the pressure and temperature dependences of the half -width 
of absorption lines are usually considered by defining a reduced path length equal to the product of 
du, the pressure, and the inverse square root of the temperature. The temperature dependence of 
the excited bands may be included by an exponential term which represents the increase of the 
population of the lower vibrational energy state with temperature (Reference 13, for example). An 
additional temperature dependence of the rotational levels, within each vibrational state, has the 
effect of altering the band shape. However, this effect is much smaller than that due to changes in 
the population of the vibrational state. The reduced path length, with these effects, is then 


du* 





ds 


(13) 


The numerical value of y was computed from Reference 14 for each spectral region, and the 
adopted values are listed in Table 1. The value /jds equals the vertical path element dh . Thus, by 
using the hydrostatic equation 

dh dP 

H “ " P ■ (14) 


the reduced path length in a direction p can be expressed by 


du*( ( ct) 




K250 " t) 


PdP 

P 2 

0 


(15) 


For the remainder of this paper, only the reduced path length in the vertical direction, u*o = l), is 
used. The reduced path length from a level P to outer space is a function of the wavelength as well 
as the temperature profile, 


U* = qH 0 



H 

0 

•v 

T(P) 


1/2 


(25O ~ T(P)) P ‘H* 


P 2 
0 


(16) 


For the numerical computations Equation 16 was reformulated in terms of the adopted coordinate 
system, 


u* 



(— - 
^250 


Vi) 


P 2 

J 


P ,’-I 


2P 0 2 


( 17 ) 


A constant temperature within each layer is assumed. The scale height H 0 for standard 
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temperature is a function of the molecular weight and, therefore, of the fractional concentration 
of each gas; adopted values are given in Table 2. * 

Since C0 2 is the primary absorbing constituent, an accurate representation of its absorption 
spectrum as a function of pressure, temperature, and path length is required. The absorption 
spectrum of C0 2 can be separated into three important regions. The near infrared region from 1 
to about 6/1 is most important in regard to the absorption of solar radiation but less significant in 
regard to the thermal emission spectrum of the clouds (Reference 15). The 9.4 and 10. 4p excited 
bands of C0 2 are within the second region considered in this calculation. They lie in a spectral 
region for which planetary observations exist (References 16-18), although they play only a minor 
role in the determination of the atmospheric structure. Any realistic interpretation of the data 
from References 16-18 must necessarily represent this part of the spectrum adequately. The 
strong bands near 15/a determine the atmospheric structure in the radiative regime; the interval of 
12-20 p is the third region considered. 

The spectrum should be divided into many intervals for accurate representation, but the number 
of intervals must be limited to avoid unreasonably long computing times. As a compromise, the 
18 intervals listed in Table 1 were selected. 


For the pressure and temperature ranges considered, the individual line shapes are controlled 
primarily by collision broadening. For high concentrations of C0 2 (q > 0.01) the distribution of 
the lines and their spacings in the major bands indicates appreciable overlapping. The square root 

law would overestimate the absorption in such 


Table 1 


Values used in the Computation of C0 2 Absorption 
for Various Spectral Intervals. 


i 

v (cm -1 ) 

m 

n 

y 

i 

0-200 

0 

0.50 

0 

2 

200-300 

0 

0.50 

0 

3 

300-400 

0 

0.50 

0 

4 

400-495 

0 

0.50 

0 

5 

495-550 

7.8 x 10 -s 

0.60 

0 

6 

550-625 

4.6 x 10 " 3 

0.42 

0 

7 

625-660 

9.0 x 10 -1 

0.38 

0 

8 

660-720 

4.9 x 10 -1 

0.42 

0 

9 

720-810 

9.5 x 10 -4 

0.40 

0 

10 

810-880 

6.5 x 10 -s 

0.55 

2200 

11 

880-920 

8.6 x 10 -6 

0.56 

2760 

12 

920-1000 

4.9 x 10" 6 

0.54 

2330 

13 

1000-1100 

4.4 x 10 -5 

0.52 

2290 

14 

1100-1400 

0 

0.50 

0 

15 

1400-2000 

0 

0.50 

0 

16 

2000-2600 

6.6 

0.50 

0 

17 

2600-8000 

1.1 x 10 -3 

0.50 

0 

18 

1700-8000 (solar) 

4.9 x 10 -7 

0.27 

0 


intervals; but the strong line approximation is 
a generalization of this law which is valid even 
when overlapping occurs. This approximation 
is expressed in terms of a single variable u* 
and is well suited to extrapolation of the ab- 
sorption to large path lengths and low pressures 
(Reference 7). 

Table 2 


Pressure Scale Height in the Atmosphere of Venus 
as a Function of the C0 2 Concentration. 


q(C0 2 per volume) 

Ho (cm) 

qH 0 (cm) 

1.0 

5.901 x 10 s 

5.901 x 10 5 

0.75 

6.491 x 10 5 

4.865 x 10 s 

0.5 

7.212 x 10 5 

3.606 x 10 s 

0.25 

8.114 x 10 s 

2.029 x 10 s 

0.1 

8.771 x 10 s 

8.771 x 10 4 

0.05 

9.015 x 10 5 

4.508 X 10 4 

0.01 

9.210 x 10 5 

9.210 x 10 3 
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For each spectral interval, the fractional absorption was computed from the expression 


A. 

1 


= 1 


(18) 


where 


T i = ( m i U *i) ni • (19) 

The m. are similar to the generalized absorption coefficients and the n. are constants determined 
by the degree of overlapping of the individual lines and the slope of the curve of growth. 

Wherever possible, the spectral intervals were chosen to coincide with the intervals used by 
Burch et al (Reference 13). The coefficients m. and n i were determined graphically by a fit of the 
strong line curve with applicable data points taken from this source. Where experimental data 
were not available the calculated tables of C0 2 absorption compiled by Stull et al. (Reference 14) 
were used. A sample is shown in Figure 2 and the various values of m and n are listed in Table 1. 

The application of the strong line approximation is valid whenever the central regions of the 
various lines become opaque. For each spectral interval a check of the validity of the strong line 



u* (cm) 

Figure 2— Illustration of strong line approximation fit for the spectral interval 720-810 cm' 1 , as derived 
from Stull et al. (Reference 14). Data points taken from Burch et al. (Reference 13) and slightly modi- 
fied are included for comparison. 
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approximation for the conditions considered was carried out. A relatively weak absorbing region 
(720-810 cm' 1 ) of the 15p C0 2 band may serve as an illustration. If the approximation holds for 
this spectral interval, then it will be valid for stronger absorbing regions. According to Plass 
(Reference 7) the strong line approximation applies for all values of /3 to an accuracy of 10 percent, 
if X > 1.63. For this interval 


/3 = 0.023 P , 

X = 0.01 y . 


( 20 ) 


According to Equations 12 and 14, u/P equals qH 0 and defines the C0 2 concentration (see Table 2). 
As can be seen, the strong line approximation is valid within an accuracy of 10 percent for C0 2 
concentrations in excess of 10' 3 and more accurate for concentrations in excess of 10" 2 , which is 
the range of interest for Venus. 


SPECIFIC INTENSITY 

The formulas derived in the section concerned with the solution of the radiative transfer 
equation must be expressed in terms of the adopted coordinate system. The specific intensity is a 
function of v, r, and /x and is conveniently expressed by upward and downward vectors. As 
mentioned before, the solution of the transfer equation appropriate here requires four components 
for the upward intensity: 

1. Thermal emission from the cloud surface 

2. Direct thermal emission from the atmosphere below interface j 

3. Downward thermal emission from the whole atmosphere, reflected upward by the clouds 

4. Downward solar radiation reflected upward by the clouds. 

The downward intensity includes values for: 

5. Direct emission from the atmosphere above interface j 

6. Direct solar radiation. 

These six components are shown in Figure 3 and discussed below. 


Surface Emission 

For a frequency interval i and an interface j , the thermal radiation emitted, which is attenu- 
ated between the surface and the interface j , is given by 


i. . 


ij 


e. B. 


(21) 
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The designation of emissivity applies directly for a solid or liquid surface, but some qualifications 
regarding the interpretation of e are necessary for a scattering medium such as a cloud surface. 
At this level the gas and cloud particles have a temperature T c . For those spectral intervals 
where the atmospheric gas is opaque, the optical properties of the cloud particles are not important, 
In regions of the spectrum where the gas between particles is transparent the optical properties of 
the particles and their concentration become very important. For these intervals c c is defined as 
a boundary value, which indicates how much smaller the emerging radiation is compared with a 
blackbody at temperature T c . The cloud emissivity was taken as unity in the strong absorption 
bands and set equal to a constant value in the transparent regions. In general, e ic may be a com- 
plicated function of wavelength and ^ but this dependence is ignored. 

Emission from the Atmosphere below Interface j 

The contribution from each layer is proportional to the Planck function, an attenuation factor, 
and the optical thickness, 
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33 


T . . 


( 22 ) 


I. . 


l J 


1 “ T. . 

! * LL 


Ar. 

i j 


Emission from the Atmosphere Reflected by the Clouds 

This component is the product of the downward flux at the cloud surface, the reflectivity, and 
the attenuation factor of the atmosphere between the surface and r . . Since constant temperature 
is assumed within a layer, the flux integral (see below) can be integrated for each layer by using 
the recursion formula for exponential integrals, 


E n (r) dr - dE n+1 (r) . 


(23) 


The reflected intensity becomes 


T ic ” T ' ' \ 1 

"" V k +1 ( T ic _ T i, /♦.) - £ 2 +1 ( T ic " T iy') 


I . - r e 

1 J 1 C 


J = 1 


(24) 


Reflected Solar Radiation 

Similarly, the reflected solar intensity is calculated from 


"18, j 


77 ^ r i8 


e U" ^ / . 


(25) 


Within the spectral range of 1.25 to 5.9/i the reflectivity of the clouds for solar radiation was taken 
equal to 0.6 which is consistent with recent observations (Reference 19). Again clouds were con- 
sidered to be perfect diffusers. 


Direct Emission from Atmospheric Layers above Interface j 

This contribution is analogous to the emission from lower layers: 

At. .* 

■ (26) 

/=i 

Direct Solar Radiation 

The direct solar radiation is given by the solar constant for 1.25 to 5.9 M appropriate for the 
distance of Venus (0.047 watt cm' 2 ): 
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(27) 


I 


18 , j 



L n 


FLUX 

The flux is derived conventionally by integration of the components of specific intensity over 
(Reference 9, for example): 


77F 


I j± dc 


2tt I I /x d/j. 


Integration leads to exponential integrals of the form of 


(28) 


1 _ 7 


/x d/x 


E 2 

— + 1 


CO 


(29) 


4 n 


7T E 2 ( t ) 


(30) 


The results are expressed by exponential integrals tabulated for integral order numbers (Refer- 
ence 20). For all wavelength intervals concerned with thermal radiation (i = 1 to 17) n is close 
to 1/2; for integration, n was taken equal to 1/2, which leads to exponential integrals of the fifth 
and fourth orders respectively. For the interval concerned with solar radiation (i = 18) the 
exponent n is 0.27; for integration, n was taken equal to 1/4, which leads in this case to expo- 
nential integrals of the ninth and eighth orders respectively. 

Analogous to the intensity, the upward flux consists of four components, and the downward flux 
consists of two. 

The individual components of the upward flux F. / are 
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4. 


(34) 




The components for the downward flux f. r are 


5 . 4 


/ ' B ii* [ E s( T ii ~ T i.i'+i) " E s( T ii ~ T ij')] 


i « l 


6- 7i‘ 


T. 


(35) 


(36) 


The net flux (-nF v = nF* - -nFj) is calculated for each layer j , for each wavelength interval i, 
and for the sum over all spectral intervals. 


CONDITION OF RADIATIVE EQUILIBRIUM 

The general condition of radiative equilibrium (Equation 8) will now be applied to an atmos- 
phere consisting of many layers of finite optical thickness. For an element of area dA in a layer 
of optical thickness A t v , the element of reduced volume is 


dV 



(37) 


In this expression f ± dA represents the projection of the element of area in a direction /z, and 
( Au */ m) is the effective element of the radiating gas. For convenience and consistency with the 
equation for r (Equation 2), i v in Equation 8 may be included with A u *. Upon inserting Equation 
37 into Equation 8, the factor (mA u *4z) n e‘ AT/M " becomes just the increment of absorptance of the 
layer. Upon integration from 0 to At the radiative equilibrium condition becomes 




_ AT X 

e ) fxdfj. dv 



e 



fj. dfi dv 


(38) 


The computed intensities were inserted into Equation 38 and integration over p was carried out 
analytically. In terms of the adopted coordinate system, the radiative equilibrium condition is 
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17 



i = 1 


The terms on the left side axe arranged in the adopted sequence; the contribution from the 
surface, then the contribution from the atmosphere (direct and reflected), and finally solar radi- 
ation, (direct and reflected). 

Radiation from the atmospheric layers above and below layer j is combined into one expression 
by the use of absolute value signs in the arguments of the exponential integrals. The right side 
represents the emission from layer j and the left side represents radiation absorbed by layer j. 

The numerical solution of Equation 39 was carried out on an IBM 7094 computer. For each 
layer j and temperature T (between 100 and 350°K) the i wave number intervals of the right side 
of Equation 39 were summed by using the appropriate value of Ar ;j . A functional relationship be- 
tweenthe total thermal emission E ; from a layer at temperature T was tabulated. The thermal 
emissionfrom a layer dependsonthe Planck function and the optical thickness; both are functions of 
the temperature. The optical thickness depends on temperature through u* (Equation 16). With 
this set of E vs. t tables for each layer, the calculation of the temperature profile continued as follows. 

First the energy absorbed by a layer j (left side of Equation 39) was computed from an esti- 
mate of the temperature profile. After the absorbed energy was obtained, the computer searched 
the previously computed table for a temperature for layer j which balanced Equation 39. In other 
words, it determined the temperature at which layer j could reradiate the absorbed energy. This 
was repeated for each layer and the resulting temperature distribution was checked for constancy 
of the net flux. If the flux criterion was not satisfied, the computed temperature distribution was 
re-inserted and the iteration procedure continued. 
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The convergence criterion used in terminating the calculation was 


r _ r 

(layer 10) (layer j) 

- 

(layer 10) 


(40) 


where a was chosen according to the desired accuracy. In general a was taken equal to 0.005. 
Obviously, higher accuracy would require more iterations. Once the convergence criterion was 
satisfied, the limb function was calculated in the final step, 


I, + (J = 1. M) 
Ii + (J = 1 . M = 1) 


(41) 


The limb function was computed for each spectral interval i and for a range of specified /j., by 
using the computed equilibrium temperature distribution. This expression (41) includes all upward 
components of the specific intensity. 

The iteration process converges rapidly for low surface pressures, where, for a particular 
layer, radiation from the surface dominates radiation from other atmospheric layers. For optically 
dense atmospheres the convergence process takes much longer. These effects are demonstrated 
in Figures 4 and 5 where samples of the iteration process for low and high surface pressures and 
for low and high initial temperature estimates are shown. More general studies of the convergence 
of iterative processes are available (Reference 21, for example). 



TEMPERATURE (°K) 

Figure 4— Sample of the convergence of the iteration process for an optically thin atmos- 
phere. The numerals adjacent to each curve in the figure indicate the iteration number 
corresponding to each temperature distribution shown. 
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TEMPERATURE (°K) 

Figure 5-Sample of the convergence of the iteration process for a more optically dense 
atmosphere than that considered in Figure 4. The numerals adjacent to each curve in 
the figure indicate the iteration number corresponding to each temperature distribution 
shown . 

RESULTS 

As stated earlier, a portion of the atmosphere of Venus was selected for analysis with the 
formulated program. Several sets of good observational material presently exist which can serve 
as a basis for comparison with the computed results. 

The atmosphere of Venus above the clouds is more amenable to analysis and consequently is 
better understood than the region between the surface and the upper cloud layer (References 22-27). 
However, starting from essentially the same set of observational data, various students of Venus 
have reached different conclusions. Estimates of the cloud top pressure range from 0.007 to 1 atm 
(References 23 and 28-30). Large discrepancies exist in estimates of the water vapor content 
(References 31-33), and the abundance of C0 2 (References 34-36) is still subject to discussion. The 
reasons for the discrepancies are many-fold, but the solution to the problem will depend on better 
observations and more refined interpretation of existing data. 

The interpretation of the intensities measured within the atmospheric window by Sinton and 
Strong (Reference 16), Murray, Wildey and Westphal (Reference 17), and by Chase, Kaplan, and 
Neugebauer (Reference 18) clearly requires some knowledge of the temperature profile in the upper 
atmosphere of Venus. So far, only the black-transparent model used by Mintz (Reference 37) and 
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the grey model used by Rasool (Reference 38) and Ostriker (Reference 39) have been applied. None 
of these calculations consider the strong temperature dependence of the 9.4 and 10.4 m C0 2 excited 
bands or the solar heating in the clear atmosphere. 

It is instructive to note that early application of the grey atmospheric model to the earth's 
atmosphere by Humphreys (Reference 40), Gold (Reference 41), and Emden (Reference 42) yielded 
only crude agreement between calculations and reality. Subsequent treatments (References 5 and 
43-47, for example) which included the wavelength dependence of the H a O, C0 2 , and 0 3 absorption 
were necessary for satisfactory temperature profiles. Similarly, a more exact and detailed study 
of the atmosphere of Venus should give better results. 

Highly developed nongrey atmospheric models (Reference 48, for example) used to describe 
the earth's stratosphere exist but have not been applied to Venus. However, recent work by Arking 
(Reference 49) and Ohring (Reference 46) with nongrey models have been applied to the atmosphere 
of Mars. For Venus, the discrepancies in estimates of composition, cloud top pressure, and other 
parameters necessary for accurate calculations have probably discouraged the application of these 
models. Furthermore, C0 2 absorption coefficients for relatively large concentrations have become 
available only recently (References 13, 14, and 50). However, it is felt that a parametric study 
with a variety of values for cloud top pressure and temperature and for C0 2 concentration should 
improve our understanding of the physical conditions in the atmosphere of Venus. Results of the 
calculations can be compared with existing and forthcoming measurements, and the range of values 
for important parameters then, hopefully, can be narrowed. 

The effect of solar heating, carbon dioxide concentration, and assumed cloud top pressure, 
temperature, and emissivity, on the atmospheric temperature above the clouds can now be dis- 
cussed. A number of conclusions can be drawn from this parametric study. However, the attempt 
to match, in detail, the calculated intensities, fluxes, and limb functions with the observations is 
still underway. 

The pronounced effect of solar heating in 
the clear atmosphere is apparent in Figure 6. 

A realistic interpretation of planetary obser- 
vations of a partially illuminated disk re- 
quires further averaging over zenith angles. 

The temperature maximum disappears for 
larger zenith angles and the temperature dis- 
tribution approaches the equilibrium condition 
on the dark side. Three factors make the 
large amount of solar heating understandable. 

Relative to the earth, the solar constant is 
nearly double; secondly, the high reflectivity 
of the clouds gives rise to strong reflected 
radiation and consequently to an additional 
heat source; and finally, the C0 2 concentration 



150 175 200 225 250 


TEMPERATURE (°K) 

Figure 6— Radiative equilibrium temperature 
distribution in the Venus atmosphere, illus- 
trating the effects of solar heating for various 
solar zenith angles. 
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is much higher. The high C0 2 concentration increases the ability of the layers to cool very ef- 
fectively, but at long wavelengths many spectral regions are already saturated. This atmospheric 
greenhouse effect is very pronounced. The near infrared absorption region considered in this 
discussion, spectral interval number 18 (1.25-5.9p), would saturate at much higher pressures 
(between 3 and 10 atmospheres). Observations at the 8-12/J. window and at 3.75 p (Reference 51) 
indicate essentially the same brightness temperatures in the sunlit and dark hemispheres. This is 
not in contradiction to the calculated large temperature differences in the clear atmosphere since 
the observations exist only for relatively transparent spectral regions. Larger differences are 
expected in moderately strong absorption regions, but unfortunately data are not available there. 

The equilibrium temperatures with and without solar heating for various cloud top pressures 
are illustrated in Figure 7. For low cloud top pressures and consequently low optical thicknesses, 
the surface can be observed more readily and hence the effective temperature approaches the cloud 
top temperature. The discontinuity in the temperature profile near the lower boundary is charac- 
teristic of radiative equilibrium considerations in a bounded, finite atmosphere and has been dis- 
cussed by several authors (References 5 and 52, for example). A more realistic representation of 
the cloud top which does not assume a "solid surface" will remove the discontinuity and limit the 
gradient to an adiabatic one. If necessary, allowance for convective as well as radiative transport 
of energy must be made. It can be also seen from Figure 7 that some of the calculated gradients 
approach and even exceed the adiabatic one near the cloud top surface. Convection would be in- 
duced and an upward transport of cloud particles would occur. 

A curve which indicates saturation temperatures for C0 2 is included in Figure 7. For the 
conditions considered, condensation and possible formation of C0 2 clouds is unlikely. The wide 
spread and the slope of temperature profiles shown in Figure 7 make it difficult to extrapolate the 
temperatures, in detail, to much lower pressures. Caution must be exercised, especially in extra- 
polations to the low pressures for the occultation measurements =2.6 X 10 " 6 atm (Reference 53). 

The effects of C0 2 concentration on the temperature distribution (Figure 8) are different for 
the dark and sunlit hemispheres of the planet. On the dark side, the temperatures depend very little 
on the C0 2 concentration, with the exception of very high layers which can cool more effectively 
for the higher concentrations. The effects of solar heating have already been discussed. 

The temperature distribution of a grey atmospheric model is included for comparison and 
clearly illustrates the difference between the isothermal character of the grey model and the 
larger gradients obtained by applying nongrey absorption. 

The final illustration (Figure 9) demonstrates the effect of nonblackness for the cloud top sur- 
face. The total net flux is kept approximately the same for all cases by adjusting the surface 
temperature upwards as the emissivity is lowered. As discussed before, the value of e was set 
equal to unity for the intervals i = 7, 8, and 16, where C0 2 absorption is strongest. In all other 
intervals £ was set equal to the value designated in the figure. The effect on the resulting temper- 
ature distribution is small. 
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Figure 9— Radiative equilibrium temperature distributions illustrating the effect of cloud top emissivity 
and temperature, for an almost constant net flux, for the dark and sunlit hemispheres. 


In summary then, radiative transfer calculations, which include pressure, temperature, and 
wavelength dependence of molecular absorption, are a strong analytical tool for the exploration of 
planetary atmospheres. The full capability of this tool has by far not been exhausted. Further 
analysis of the spectral and angular dependence of the emitted energy is required for more specific 
conclusions. 
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